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We discuss the behavior of the entanglement entropy of the ground state in various collective sys- 
tems. Results for general quadratic two-mode boson models are given, yielding the relation between 
quantum phase transitions of the system (signaled by a divergence of the entanglement entropy) and 
the excitation energies. Such systems naturally arise when expanding collective spin Hamiltonians at 
leading order via the Holstein-Primakoff mapping. In a second step, we analyze several such models 
(the Dicke model, the two-level BCS model, the Lieb-Mattis model and the Lipkin-Meshkov-Glick 
model) and investigate the properties of the entanglement entropy in the whole parameter range. 
We show that when the system contains gapless excitations the entanglement entropy of the ground 
state diverges with increasing system size. We derive and classify the scaling behaviors that can be 
met. 
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I. INTRODUCTION 

One of the most striking aspects of quantum mechan- 
ics is certainly the superposition principle which has no 
counterpart in classical physics. A direct consequence 
of this principle is the entanglement of physical states 
whose characterization has been a subject of intense re- 
search in the recent years. However, the interest for the 
entropy in quantum systems has first emerged in black 
holes physics p[ pre-empting the recent favour in con- 
densed matter physics. The relationship between entan- 
glement and quantum phase transitions (QPTs) has es- 
pecially drawn much attention since the original works 
on one-dimensional (ID) spin systems [3, U0, @- ^ n 
these studies, it has been shown that entanglement mea- 
sures are very sensitive to the presence of a critical point 
in the phase diagram. For example, in the ID quantum 
Ising model in a transverse field, the entanglement en- 
tropy of the ground state has been shown to be finite 
everywhere except at the critical field where it diverges. 
This entropy, which is the central topic of this paper is 
defined as follows. Suppose that we split the degrees of 
freedom of the system under consideration into two sub- 
systems A and B. The entanglement entropy £ of any 
pure quantum state p — \il))(ip\ with respect to this bi- 
partition is defined as the von Neumann entropy of the 
reduced density matrix 

£ = -Tta (pa In p A ) = -Tr B (p B lnpg) , (1) 

where pa,b = ^B,A P- Several fundamental questions 
arise concerning this entropy. In particular when, why 
and how does it diverge? In d-dimensional non-critical 
(gapped) systems, the entanglement entropy is generi- 
cally proportional to the surface area L^ 1 of the consid- 
ered subsystem with linear size L (area law). In critical 
systems, logarithmic corrections to this scaling may oc- 
cur. In particular, a logarithmic scaling £ oc InL is found 
in critical one-dimensional boson and fermion (and thus 



also spin) systems d S S S B i, [M HH, El ■ Note that 
in conformal field theory, the prefactor is determined by 
the central charges fl3l [l4| . In higher-dimensional crit- 
ical systems, the scaling behavior is different for bosons 
and fermions. For bosons, field theoretical and numerical 
results show that the area law prevails also in the critical 
case p], [H, [H, [I?], EH- For fermions, the area law is in 
this case corrected by a logarithmic factor pl| 
Note also that recently, noncritical fermionic systems 
with rapidly decaying interactions have been shown to 
violate this area law [2l| . 

All the above cited studies focus on low-dimensional 
systems. The aim of the present work is to investi- 
gate collective systems which, in a sense, can be seen 
as infinite-dimensional. The physics of these models can 
be easily analyzed by a classical approach but the en- 
tanglement entropy comes from the quantum fluctua- 
tions around the classical ground state. Here, we map 
those collective systems onto models with two interact- 
ing bosonic modes, each of them corresponding to one 
part of a bipartite splitting of the original system. We 
show that the scaling of the entanglement entropy de- 
pends mainly on two factors: the number of Goldstone 
modes (0, 1 or 2) or, when approaching a critical point, 
the number of vanishing excitation gaps (1 or 2). After 
giving the general results for quadratic two-mode boson 
Hamiltonians, we discuss in detail several collective mod- 
els covering different behavior of entropy scaling. 

The structure of this paper is the following. In Sec. 
im we derive the expression of the ground state entan- 
glement entropy for general quadratic two-mode boson 
Hamiltonians and extract its explicit dependence on the 
excitation energies. Section Hill is devoted to the study of 
the Dicke model for which the entropy is shown to diverge 
as i In N at the critical point where N is the number of 
two-level systems. In Sec. IIV1 we analyze the two-level 
BCS model in which the entropy diverges as | In N in 
the broken phase and as |ln./V at the critical point. In 
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the antiferromagnetic long-range Lieb-Mattis model dis- 
cussed in Sec. [V] we exactly compute the ground state 
entanglement entropy which behaves as In TV. Section 
IVII focuses on the Lipkin-Meshkov-Glick (LMG) model 
for which we analyze the dependence of the entropy on 
several parameters (anisotropy, system and subsystem 
sizes). In each case, we compare our analytical predic- 
tions based on the scaling hypothesis (see Appendix [AJ 
with numerics. 



particle energies are necessarily positive. This is given, 
if we restrict to u)± > 0. After a Bogoliubov transforma- 
tion, H is diagonalized in terms of new bosonic operators 
b O, 



H = (P) T A 6 + i Tr(A — V) , 



and the Green's function matrices read 



(5) 



II. GENERAL QUADRATIC TWO-MODE 
BOSON MODEL 



G++ = i? t (^)w 1/2 i? t (V)A~ 1 i?W^ 1/2 i?(^), (6a) 
G = -R\<p) uj- 1/2 R^ip) A Rty) lu- 1 / 2 R(ip) . (6b) 



Let us consider the most general quadratic two-mode 
boson Hamiltonian 



H = (tf) T Vb + ^ (b T )' J Wb 1 + h 



l r 



(2) 



where b = (6i, 62)- We restrict here to the case where the 
coupling matrices V and W are real so that there are six 
free parameters. The aim is to find a simple expression 
for the entanglement entropy between both modes. As 
shown in Ref. 18], the reduced density matrix obtained 
by tracing over one of the two modes can be expressed in 
terms of the Green function matrices defined as G^ + = 

((b\ + h){b] + b-)\ and Grr = ((&]_- &,.)(&} - bj)). _ It 
is therefore convenient to reparametrizc the Hamiltonian 
© as 

V -W = R\ip)wR{ip), (3a) 
V + W = R\ip) or* R\4>) A 2 R(ip) uj-^R((p) , (3b) 

where 



R{<p) 



cos <p sm ip 
— sin ip cos ip 



(4) 



u; = diag(w + ,w_) and A = diag(A + , A_). The matrices 
V ± W need to be positive definite, because the single- 

I 



The entanglement entropy between modes 1 and 2 is then 
given by [181 ] 

M+l m+1 m-1 

2 2 2 2 ' V ' 

l l = \f -GfiG 11 = y / -GjjG 22 . (7b) 

Now, let us express all quantities in terms of the initial 
coupling matrices. Therefore, let us set: 



V — W 



do d 
d di 



which leads to 



d Q + d x ± eV4d 2 + (do - di) 2 



where e = sign(d — di)- 
Similarly, setting: 



s s 

s Si 



V + W = 
one obtains the excitation energies: 



(8) 



(9) 



(10) 



with 



A H 



1 r 



2ds + c?qSq + d\S\ ± e'-\/ (2ds + doso + disi) 2 — 4(do<ii — d 2 ){soS\ — s 2 ) 



1 _ . f e2d[s(dp + di) + d{s + si)] + (d - di){d s - diSi) 
" " 8181 V 4d2 + (do - d,) 2 



(11) 



(12) 



Of course, A± have to be positive, giving the parameter range for which this transformation is well-defined. It is also 
useful to introduce the two parameters 



t = taa(2tp) = 2d/{d - dx) 



(13) 



u = tan(2^) 



2{dpd 1 - d^yfjsjdo - dx) - d(s - Sl )} 
2d[s(d + dx) + d(s + Si)] + (d - dx){d so - dxsx) 



(14) 
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The quantity /i denned in Eq. (]Tb[) is then given by 



(15) 



with 



_ t(l + M 2 ) 1 / 2 (^-Q±[2 M + ^ + ^ 1 )] 



4(l+i 2 ) 1 /2(l +u 2)l/2 



and 




(17) 



Note that these expressions are still valid for the special 
case do = di provided one sets e = sign(d). 

This simple form allows one to investigate the various 
behaviors of the entropy that can be met in two-mode 
boson systems. Indeed, Eq. (TT5"]) clearly states that the 
entropy diverges when any of the quantities ui±, A± van- 
ishes or diverges and several cases must then be distin- 
guished. 

In the following, we will only focus on systems for 
which gaps vanish, signalling the presence of a quantum 
phase transition. In the next section, we consider the 
Dicke model [23| in which one gap vanishes at the crit- 
ical point but the other remains finite. In Sec. IIV1 we 
consider the two-level BCS model for which both gaps 
vanish simultaneously at the transition. We show how 
these different behaviors lead to different finite-size scal- 
ing for the entropy at the critical point. 



III. THE DICKE MODEL 

Let us consider the single-mode Dicke Hamiltonian [23[ 
which describes the interaction of N two-level systems 
with a single bosonic mode 

H = Lu S z +Lu^a + ^(a^ +a){S+ + S-) , (18) 



where and a are bosonic creation and annihilation op- 
erators satifying [a,a^] = 1. The collective spin operators 
are defined as S a — J2iLi a a/^ where the cr Q 's are the 
Pauli matrices, and S± — S x ± i S v . 

We refer the reader to Ref. [24[ for a description of 
the phase diagram as well as the super-radiance physics 
in this model. Here, we focus on the entanglement en- 
tropy studied in Refs. 0, [26[ and discuss the linear en- 
tropy recently analyzed in Ref. [13] • Let us also mention 
that other entanglement measures such as the concur- 
rence have been studied in this model [25|, [H, [29| . 

Our goal here is not to reproduce the results of Lam- 
bert et al. [25| but to show how to extract the finite-size 
behavior of the entropy at the critical point. Thus, we 
only consider the normal (symmetric) phase (A < A c = 
^/ojujq/2) characterized by \impf^ 00 (a'a) /N = 0. The 



Holstein-Primakoff (HP) boson representation of the an- 
gular momentum [30j reads: 



tfN l ' 2 {l - tfb/N) 1 ' 2 = {S-) f 



(19) 
(20) 



with [b, &t] = 1 5 so that we now have to deal with a two- 
boson problem. In the thermodynamical limit and in the 
normal phase, one has (b'b)/N -C 1 and we can expand 
the square root in Q19p to obtain the following form of 
the Hamiltonian: 



H = -y^o + uotfb + Loa)a + A (a) + a) (b*< + b) 



-0{l/N) 



(21) 



At this order, one thus has to deal with a quadratic 
Hamiltonian and we can make use of the results given in 
Sec. |TTJ As found in Ref. [26j, the two excitation energies 
are given by 



\Ju 2 + us 2 ± yT~6A 2 wwo + (w 2 — uJqY 

71 



.(22) 



Note that A±, as defined in Eq. (jlip. depends on 
e = sign(w 2 — Wq) which is assumed here to be pos- 
itive without loss of generality. At the critical point 
A c = t/lJujo/2, the system undergoes a second-order QPT 
and A_ vanishes as (A c — A) 1 / 2 whereas A + = \/ui 2 + u)q 
remains finite. 

In this model, the parameter fi is given by: 



,(0) _ 



AloloqX 2 



16A 2 wu;o + (uj 2 — u>q) 



A_ 
A" 



- 2 



where the superscript " (0)" refers to the order in the 1/N 
expansion. The entanglement entropy between modes a 
and b is then given by Eq. (|7a|) . In the vicinity of the 
critical point, the entropy thus diverges as: 



£(°)=l+ln 



A 



2(A C -A) (oj 2 +lo 2 Y 



1/4 



+ 0[(Ac-A) 1/2 ]. 

We note that this expression is in agreement with those 
obtained in Ref. [26| with a different approach. However, 
the constant term in this latter study is wrong and we 
give here the correct value. 

To extract the behavior of the entropy at the critical 
point, one simply invokes the scaling hypothesis discussed 
in Ref. [29] for this model and detailed in Appendix [Al 
In the present context, one indeed has, near the critical 

point, jln(A c - A) which yields £ - \\nN. 

Note that this exponent 1/6 is close to the numerical 
result (0.14 ±0.01) obtained from exact diagonalizations 
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Similarly, one can easily compute the linear entropy 
which reads : 



c(0) _ , _ 1 
C lin. — 1 



(0) ' 



(24) 



as shown in Appendix [Bj In the vicinity of A c , one has 



AO) _ j _ 



(A c -A) {oj 2 +ujIY 
8A5 



1/4 



+ C[(A C -A) 3 / 4 ] 



Consequently, using the scaling hypothesis, and 
AO) 



"Tin. 



(25) 
since 

1 - (A c - A) 1 / 4 one expects £ Un . - 1 ~ A" 1 / 6 at 
the transition as predicted in Ref. [27| in the adiabatic 
regime (u>/loq <C 1). Here again, we wish to emphasize 
that the expression given by Liberti et al. is only valid in 
this regime. Thus, their expression of the linear entropy 
is only correct at first order in lu/wq as can be easily 
checked by expanding (|25[) in powers of this parameter 
uj/ujq. 

To conclude this section, we mention that the same 
calculation can be performed in the broken phase A > A c 
(see Ref. [26]). Further, the exact correction of order 
1/N could, in principle, be easily calculated but its ex- 
pression would be too long to be given for this model. We 
shall compute it for the systems considered thereafter. 



IV. THE TWO-LEVEL BCS MODEL 

We now turn to another problem which can also be 
mapped onto a two-mode boson system, but where both 
gaps vanish at the critical point. The two- level BCS 
model originates from superconductivity [3lL [32j but can 
be translated into a simple spin model (see for instance 
Ref. 33]). Its entanglement properties have recently 
been discussed but only from the concurrence point of 
view pjil [35| . Here, we shall show that its entropy cap- 
tures some nontrivial features due to the vanishing gaps. 

The two-level BCS Hamiltonian reads 



H 



Sf, 



KS?-S?) 



(26) 



Here, the indices A and B refer to two distinct sets of 
A/2 spins 1/2 subjected to antiparallel magnetic fields. 
Here, we have introduced the collective spin operators 
S£' B = YlieA,B °a/2j and S a = S£ + S% (where a = 
x,y,z). We further restrict our study to the total mag- 
netization S z = which is reminiscent of the fermion 
number conservation in the initial BCS problem. Let us 
first discuss the symmetric phase h > 1 characterized by 
liniAr^ 4(Sf)/N = - limbec 4(Sf )/N = 1. As previ- 
ously, we perform a 1/N expansion using the HP boson 
representation of the spin operators: 



s* 



N/4- a f a, 



tfb - A/4, 



(27) 



{S A y = (A/2) 1 / 2 (1 - 2a t a/A) 1/2 a , (28) 



(S*)* = (A/2) 1 / 2 6 t (1 - 2tfb/N) 



1/2 



(29) 



with Si' B = S^ B ±iSf' B and [a, at] = [6, &t] = 1. I n the 
thermodynamical limit and in the normal phase, one has 
(a^a)/N < 1 and {tfb)/N < 1 so that the Hamiltonian 
can be written as: 

H = -^-i+(fc-i)( to + 6t 6 )_l( t 6 t + o6) 

+ — ( a) a + b^b + a? a + b^b 2 
N V 

+ a) a b + a 1 a6M 
2A V / 

+C(1/A 2 ). (30) 

At order A , the Hamiltonian is easily diagonalized and 
one then finds that, in this phase, both gaps are equal to 
A = \Jh(h — 1). Thus, at the critical point h = 1, both 
gaps vanish (k = 1). 

The parameter fi can however still be expressed in 
terms of the gap 



(o) _ 



2/i-l 
2A ' 



(31) 



where the superscript " (0)" again refers to the order in 
the 1/N expansion. The entanglement entropy between 
sets A and B, in the thermodynamical limit and in the 
symmetric phase is then given by Eq. (|7a|) . In the vicin- 
ity of the critical point, the entropy thus behaves as: 



£ (0) = 1 - 21n2 -ln(h- 1) 1/2 + 0(h - 1) 



(32) 



The scaling hypothesis (see Appendix [XJ immediately 
implies £ ~ | In A at the transition point. 

Before comparing with numerics and discussing the 
broken phase, we also give the next order correction 
which can be computed exactly in this model. At or- 
der 1/N, the Hamiltonian can be diagonalized using the 
canonical transformations method [291 ]. This allows one 
to evaluate the correction of order 1/N to /j, 



(i) 



h 1 

a7 + a7 



(33) 



As explained in Ref. (36|, p/ 1 ^ is sufficient to determine 
the first correction to the entropy (see Appendix The 
entropy at order 1/N is indeed obtained by Taylor ex- 
panding expression (|7a[) where /x is also computed at this 
order as above. One then obtains: 



*(°) + 1 



(o) 



(34) 



This expression allows to check explicitly the scaling hy- 
pothesis (see Appendix E|) since, at this order, and in the 
vicinity of the critical point, one has 



e = e +N 



(h-iy 1 ' 2 



1 - 



N(h-1) 3 / 2 



(35) 
(36) 
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FIG. 1: (Color online) Entanglement entropy in the two-level 
BCS model as a function of h for N = 32, 64, 128, 256 (from 
numerics [black lines]) and oo (£^°' [red line]). Arrows indi- 
cate the behavior of the finite-size correction in various re- 
gions. 





FIG. 3: (Color online) Entanglement entropy as a function of 
In N in the two-level BCS model for h = 0, 1 /2, 1 and analytic 
predictions: slope 1/2 (green line) for h < 1, slope 1/3 (blue 
line) for h = 1. 



whose entropy is known [3?], HI] to diverge as |lnJV. 
Since no drastic physical changes occur when varying the 
field below ft, < 1, it is reasonable to expect a similar 
behavior in the entire broken phase. We have checked 
this handwaving argument with numerics and, indeed, 
the ground state entropy diverges as i In N in the entire 
broken phase (see Fig. [3]). 

To complete this study, it is interesting to analyze a 
similar model for which the classical analysis predicts 
two-Goldstone modes instead of one. The Lieb-Mattis 
Hamiltonian 139] is a natural candidate. 



FIG. 2: (Color online) Behavior of AT[£ num — £ (0 ^] as a func- 
tion of h in the two-level BCS model for N = 32, 64, 128, 256 
(from numerics [black lines]) and oo (£^ [red line]). 



As can be seen in Fig. [T] an excellent agreement is found 
between the entanglement entropy for increasing N and 
the thermodynamic limit. The leading finite N correc- 
tion is also extracted from numerics and compared to the 
analytical result (see Fig. [2]) given above confirming the 
validity of the scaling hypothesis. 

Let us now turn to the broken phase h < 1. In this 
phase, at order N°, one of the two gaps vanishes whereas 
the other remains finite. More precisely, one has 



A + = y/1 - V . 

A_ = 0. 



(37) 
(38) 



We do not give here the details of this analysis since it 
is similar to the isotropic LMG model discussed in the 
next section. Unfortunately, one can not use the scaling 
hypothesis to determine the behavior of the entropy due 
to the presence of the Goldstone mode in this parameter 
range (A_ = 0). However, one can exactly determine the 
entropy in zero field. There, the ground state is given by 
the Dicke state corresponding to S = N/2 and S z — 



THE LIEB-MATTIS MODEL 



The Lieb-Mattis Hamiltonian [39j is given by 



H = - S A ■ S B = — 

N 2N 



S 2 - S A Y - s l 



(39) 



As for the two-level BCS model, the indices A and B 
refer to two distinct sets of N/2 spins 1/2. The (non- 
degenerate) ground state |^o) is readily found to be 

the eigenstate of |s 2 , (S- 4 ) 2 , (S 8 ) 2 j, with eigenvalues 



{0, f (f + 1), f (f + 1) }• However, the classical ground 
state is found to have two Goldstone modes stemming 
from the isotropy of the Hamiltonian (|39]l. This suggests 
that the scaling behavior of the entanglement entropy 
here should differ from those of the models considered 
previously. As in the broken phase of the BCS model, the 
diagonalization of the quadratic Hamiltonian that would 
derive from a HP expansion does not give the correct 
ground state because of the Goldstone modes. Neverthe- 
less, the entanglement entropy of the ground state can 
be evaluated exactly since: 



1 



N/4. 

= E (-D 

1 M=~N/4 



-M+iV/4 



\M,-M), (40) 
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where the state \M,—M) denotes the eigenstate of 
{(S A ) 2 ,S A , (S B ) 2 ,S? } with eigenvalues {f(f + l),M, 

^ (-j + 1) , M | . The reduced density matrix pj( obtained 
by tracing over set B is thus easily computed and reads 
in the eigenbasis of j^S* 4 ) 2 , S^j 



PA = Tr B \^o)(ip \ 



N 



-\M)(M\. (41) 



The entanglement entropy is thus given by 

£ = -Trplnp = ln(f + 1), (42) 

and thus behaves as IniV at large N. 

To conclude this brief review of entanglement entropy 
in collective systems, we shall now consider a collective 
spin model with an arbitrary bipartition. 



VI. THE LIPKIN-MESHKOV-GLICK MODEL 



where S a = '}2 i & l a /2, an d the a a 's are the Pauli matri- 
ces. It describes a system of N spin 1/2 mutually inter- 
acting in a field h, transverse to the interaction directions 
(XY). Here, we restrict our study to the most interest- 
ing ferromagnetic case and, without loss of generality, we 
assume ^ 7 < 1 and h ^ 0. The isotropic case 7 = 1 is 
discussed separately at the end of the section. 

As was early discussed in the original work of Lipkin, 
Meshkov and Glick, this system undergoes a second-order 
QPT at h — 1, between a symmetric (h > 1) and a 
broken (h < 1) phase. In both regimes, the ground state 
lies in the maximum spin sector S — N/2. To analyze the 
entanglement entropy, we partition the set of N spins into 
two blocks A and B, of sizes L and (N — L) respectively 
[FH ] . We also introduce the corresponding spin operators 

Sot' —^ieA.B^a/^- 

Let us first focus on the symmetric phase h > 1 in 
which the classical ground state is fully polarized in the 
z-direction. As previously, it is convenient to use the HP 
representation of the spin operators 



Initially proposed fourty years ago by Lipkin, Meshkov 
and Glick [i^, El, |42j to describe phase transitions in nu- 
clei, the LMG model has, since then, been widely used 
to describe many physical systems such as Bose-Einstein 
condensates [43(, Josephson junctions or long-range in- 
teracting spin systems [3 [45| . Here, we shall consider 
the latter formulation which is especially well-suited to 
our problematics. The entanglement properties of the 
LMG model have drawn much attention in the last years 
HI, H S3, H H HI but the entropy has been analyt- 
ically investigated only very recently [5l| following the 
numerical work of Latorre et al. [38] . 

The Hamiltonian of the LMG model is given by 



H = - — (S 2 X 



hS z 



(43) 



S A 
S B 



L/2-a)a , 

L 1 ' 2 at (1 - cJa/L) 1 ' 2 = (S A )^ , 

(N-L)/2- 6+6, 

(N - L) 1 ' 2 (1 - 6+6/(JV - Lf/ 2 



(44) 
(45) 
(46) 
(S*)+ , (47) 



with 5± ,B = S A ' B ±iS A ' B . This transformation maps the 
LMG Hamiltonian onto a system of two bosonic modes. 
At fixed t = L/N, the Hamiltonian can be expanded in 
powers of 1/N since, in this phase, one has (a'a)/L <C 1 
and (tfb)/(N - L) < 1. At order (1/N) and for h > 1, 
one then gets 



H 



Nh I + 7 2/1-7 - 1 f + 



+ ' Zn 7 i fcta + fttfc) + 2_l[ T ( a f 2 + a 2 ) +(l- r )(V 2 + fr 2 ) +2vV(l-T)( a t6t +afe ) 



1-7 



(1 - t) (a)ab + a^ab^ + r (atfb 2 + a+6+ 2 &)] + (a+a + tfb + a+V + b^b 2 + 20+06+6) + 



1-7 



a) 2 + a 2 + 6 f2 + b 2 + 2( a^a 3 + a f3 a + b ] b 3 + 6+ 3 6 



0(1/N 2 ). 



(48) 



At order N°, this Hamiltonian is easily diagonalized and 
the excitation energies are given by: 

2/1-7-1 



(49) 
(50) 



= y/(h-l)(h-7). 

Thus, at the critical point, A_ vanishes whereas A + re- 



mains finite except for 7 = 1 (see discussion at the end 
of the section). The parameter \x is then given by: 



(o) _ 



l + r(l-r) 



h-1 

h — 7 



1/4 



h — 7 
h-1 



1/4- 



(51) 
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FIG. 4: (Color online) Entropy as a function of h in the LMG 
model for 7 = 1/4, r = 1/4 and N = 32,64,128,256 (from 
numerics [black lines]) and 00 (S^ red line). Arrows indicate 
the behavior of the finite-size correction in various regions. 
For h < 1, we plotted £ (0 ^ + In 2 to take into account the 
parity-broken nature of the classical states taken as starting 
point of the 1/TV expansion [26l . [5lJ. The inset is a zoom 
around h = ^7 where £^ ' — 0, for TV = 64 (black line) and 
00 (red line). 



which is equivalent to the form presented in Ref. [Hlj . 
The entropy is then given by Eq. (|7a|) . Figure H] shows 
a comparison between the entropy at finite TV obtained 
from numerical diagonalization of the Hamiltonian and 
the entropy £^ obtained from Eqs. (|7a|) and (f5Tj) . As 
can be seen, the entropy is finite for all h 7^ 1 and diverges 
at the transition point. 

Near the critical point h = 1 + , it behaves as 



£(0) 



~ ln(ft - 1) + i ln[r(l - r)] + \ ln(l - 7) + 



1 -ln2 + C[(/i- 1) 1/2 ] 



(52) 



where the superscript " (0)" refers to the order in the 1/TV 
expansion. As in the Dicke model studied in Sec. IIII1 the 



scaling hypothesis (see Appendix yields £ ~ g In TV. 
As can be seen in Eq. (j52|) . the entropy also diverges 
when 7 goes to 1 and when r reaches its extremal val- 
ues or 1. This suggests a possible dependence of the 
scaling variable on these parameters. To obtain this de- 
pendence, one needs, at least, to compute the correction 
of order 1/TV to the entropy. As explained in Ref. [36j, it 
is obtained by Taylor expanding expression (|7aj) at order 
1/TV where /x is also computed at this order (see Ap- 
pendix [C]). This latter step requires to diagonalizc H at 
order 1/TV. Usin g th e canonical tranformation method 
described in Ref. [29| , one then obtains the correction of 
order 1/TV to /z 



(1) _ r(l-r)(l- 7 ) 



- r 47-/i(l + 2/i + 7) h 
2A 7 + 



(53) 

As already stated in Sec. IIV1 the correction of order 1/TV 




FIG. 5: (Color online) Behavior of TV[£ num — as a func- 
tion of h in the LMG model for TV = 32, 64, 128, 256 (from 
numerics [black lines]) and 00 (S^- 1 ' red line). 



to the entropy then reads 



,(!) 



£(i) = ^i n 



(0) + i 



(o) _ 1 



(54) 



We emphasize that this is the exact correction of order 
1/TV whereas in Ref. [5l[, onr y its expansion near h = 1 
was given. We display in Fig. [5l a comparison between 
the numerical correction of order 1 /TV and the analytical 
expression obtained from Eqs. ([55]) - ([54")) . 

In the vicinity of the critical point h = 1, one thus has, 
at this order: 



c 



(i) 



(55) 



[r(l-r)] 



1/2 



1/4 



(1-7) 



1 



3(1-7) 1/2 



8N{h- I) 3 / 2 



This suggests that the scaling variable is indeed a func- 
tion of 7 but has no dependence on t. Following Ref. 
[5ll ] , we reformulate the scaling hypothesis for this model 
by assuming that in the vicinity of the critical point, a 
physical observable $ can be written as the sum of a 
regular and a singular contribution, 



where 



TV"* 



(56) 



(57) 



and £ = N(h — 1) 3 / 2 (1 — 7) -1 ^ 2 is the scaling variable. 
The exponents £j and n$ are characteristics of the 
observables $. Thus, at the critical point, one expects 

a»(0 ~ r 2 ^ /3 , so that 

= 1, 7) ~ TV-("*+ 2 ^/ 3 ) (1 - 7)«2+«i/3 _ (5g) 

We emphasize that the scaling form (|A3|) . can be 
checked up to high order from the expansion of the 
ground state energy, the gap, the magnetization, and the 
spin-spin correlation functions given in Ref. [501 ]. 
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FIG. 6: (Color online) Exponents Xi an d Xt as a function 
of 1/JV obtained from numerical diagonalization of H. For 
clarity, we plotted 2\t and 6X7 which are expected to be equal 
to 1 in the thermodynamic limit (dotted lines are guides for 
the eyes). 
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FIG. 7: (Color online) Entropy as a function of In N for 7 = 
and t — 1/2. The blue line has a slope xn — 1/6. 



The scaling law ([58]l applied to the observable e £ leads 
to the following critical behavior of the entropy 

£(/i = 1)~Xw1hJV+x 7 1h(1-7)+XtMt(1-t)], (59) 

with xn — X-y — 1/6 and Xr = 1/2- As can be seen 
in Figs. [5] and an excellent agreement between these 
predictions and numerics is observed when increasing the 
system sizes. 

Such a behavior has striking similarities with the ID 
XY model @, S 0, 1, El HI . This is due to the fact that 
the ground state of the LMG model lies in the maximum 
spin sector S — N/2. As a consequence, the relevant 
subsystem Hilbert space is spanned by the L + 1 eigen- 
states of {(S A f, Sf} with S A = L/2. The entanglement 
entropy as a function of the subsystem size L can thus 
not exceed the maximum value ln(i + 1). However, note 
that the physics of the LMG model is very different from 
ID models. 

This surprising analogy with ID XY model carries over 
to the broken phase. Indeed, at order N°, one can com- 
pute the entropy in the same way as in the broken phase. 



Therefore, one first has to perform a rotation around 
the y-axis to bring the z-axis along the classical ground 
state magnetization S cL = N/2(±sm9o,0,cos9o) with 
9q = arccos h [5Cj • To use the HP mapping, we thus need 
to rotate the spin operators according to 




cos #0 sin 6*o 

1 
— sin 8q cos 9q 




(60) 



Note that this choice corresponds to study quantum fluc- 
tuations around one of the two classical minima (see 
for discussion). Next, one uses the HP representation of 
S which is exactly the same as in Eas.(|35 ]) -(|17 |) . The 
Hamiltonian at order N° then reads: 



H 



N{l + h 2 ) h 2 + 1 



7 



7' 



4 

2V7(T 



6 f +ab 



• il r,[h +b ) r 

-0(1/N). 



(61) 



As can be easily checked, this Hamiltonian coincides with 
(|48|) for h = 1 as it should for a second order QPT. Using 
the general results of Sec. HIl one straightforwardly gets 
the two excitation energies 



A+ = 



2 - h 2 - 7 



/i 2 )(l- 7 ), 



(62) 
(63) 



and 



„<"> = 



\ 



l + r(l - 



1-h 2 
1-7 



1/4 



7 



1/4 



(64) 

It is worth noting that for h — ^7, one has = 1 so 
that the entropy vanishes (see Fig. 0]). At this point, 
the ground state (which is two-fold degenerate in this 
phase) can indeed be chosen as separable I50| . Such a 
point is often present in anisotropic systems [53j. 



The isotropic case 

For 7 = 1, the above expression of the gaps are still 
valid in the thermodynamical limit and one recovers 
the Goldstone mode predicted by the classical approach. 
This gapless excitation prevents from using the scaling 
hypothesis to determine the behavior of the entropy. For- 
tunately, the isotropic case is trivial since the eigenstates 
of the LMG Hamiltonian (|43[) are then the eigenstates 
\S,M) of S 2 and S z . The ground state is obtained for 
S = N/2 and M < N/2 for h < 1, whereas for h > 1 
one has M = N/2. Thus, in the symmetric phase, the 
entropy vanishes since the state \N/2, N/2) is separable, 
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but in the broken phase, the entropy of the Dicke state 
\N/2,M < N/2) is known to diverge as |lniV [HHl]. 
This means that the entropy is discontinuous at the crit- 
ical point h = 1 (as the concurrence [5(|) although this 
transition is often considered as second-order. 



VII. CONCLUSION AND PERSPECTIVES 

The results presented in this paper show that the en- 
tanglement entropy in collective models displays a large 
variety of behaviors. First of all one has to distinguish 
with respect to the number of Goldstone modes in the 
system (0, 1 or 2). 

1. The classical ground state has a finite degeneracy, 
i.e. no Goldstone mode: this situation is met in the 
broken phases of the Dicke and anisotropic LMG 
model, as well as in the symmetric phases of the 
models studied here. In this case, and away from 
critical points, both excitations around any classi- 
cal states are gapful and the entropy is finite (pos- 
sibly vanishing). 

2. The classical ground state has one Goldstone mode: 
this situation is met in the broken phases of the 
two-level BCS model and the isotropic LMG model. 
In this case, one excitation energy is finite and the 
other vanishes and the entropy diverges as h In N. 

3. The classical ground state has two Goldstone 
modes: this situation is met in the antiferromag- 
netic Lieb-Mattis model. In this case, the entropy 
diverges as In AT. 

In case 1, where we have no Goldstone mode in the sys- 
tem, the collective spins, S A and S B , of the two sub- 
systems may be expanded in 1/N around their classical 
orientation via Holstein-Primakoff transformation. The 
system is thus mapped to a two-mode boson model. The 
general analysis of the order iV° (at which the effective 
bosonic models are quadratic) in Sec. [Til shows how the 
entanglement entropy diverges when excitation gaps van- 
ish. Its scaling can thus be classified according to the 
number vanishing of gaps, when a critical point is ap- 
proached. In the anisotropic LMG model as in the Dicke 
model, there is only one vanishing gap and the entropy 
diverges as | In N in both models. 

In the two-level BCS model, there are two vanishing 
gaps and we obtained an entropy that diverges as ^ In N. 
However, these behaviors are found using the scaling hy- 
pothesis which relies on the scaling variable N(g — g c ) a 
(see Appendix |A|) . In models studied in this paper, one 
always has a = 3/2 due to their relationship with a quar- 
tic oscillator (54|. Nevertheless, this value depends on the 
Hamiltonian and would be different if, for instance, at the 
critical point, the terms of order A^ and 1/N vanish si- 
multaneously. In this case, the Hamiltonian would be 
sextic instead of quartic in the creation and annihilation 



operators and the scaling of the entropy would be differ- 
ent [55[ . Note also that in the isotropic LMG model, the 
entropy vanishes at the critical point even though there 
are two vanishing gaps at the critical point. 

The analysis performed here for simple models corrob- 
orates the results obtained in many ID systems where it 
has been also found that diverging entropies are always 
associated to gapless excitations |4|, |5|, [D, |8|, |l3j, LL21 . This 
should constitute a strong motivation to further inves- 
tigate collective models which are simpler to treat than 
low-dimensional ones and already exhibit nontrivial be- 
haviors of the entanglement entropy. 
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APPENDIX A: THE SCALING HYPOTHESIS 

This appendix is dedicated to a basic description of 
the scaling hypothesis. This hypothesis is useful in col- 
lective models where a 1 /N expansion is meaningful such 
as those considered in this paper. It relies on the fact 
that, in these models, the 1/AT expansion of the expec- 
tation value of any physical observable can be written as 
the sum of a regular and a singular contribution, 

^ N (h n ) = ^(g) + ^(g). (Al) 

Here, singular means that the function and/or its deriva- 
tives with respect to the control parameter g driving the 
transition diverge at the critical point g Cl following a 
power law. In addition, one has 

*7 S (9 ~ 9c) ~ (g ~ffV g[Ar(g - g D )°] . (A2) 

The exponents f$ and n$ are characteristics of the ob- 
servables $ but the exponent a depends on the model 
considered as early understood by Botet and Jullien [45[ . 
For the models considered in this paper, we have explic- 
itly checked that a = 3/2 which is reminiscent of the 
behavior of the gap in the vicinity of the critical point 
which vanishes as (g — gc) 1 ^ 2 ■ 

The behavior of <! > ™ s at the critical point g c is then 
obtained by noting that since no divergence can occur 
at finite N, one must have !F$(x),~ x ~^/ a and conse- 
quently 

^% ng {g = g c ) ~ AT-(«*+«Wa) _ ( A3 ) 

Finally, since we are mainly interested here in the en- 
tropy, let us simply note that to obtain a power-law 
behavior, the observable to be considered is e £ which 
indeed bevaves as suggested in (|A2[) . As a result, if 
£ ~ — £,£ hi(<7 — g c ) near the critical point, the scaling 
hypothesis predicts that it diverges as — In A". 
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This scaling hypothesis has been checked up to high 
order in the several models and allowed us to obtain the 
finite-size scaling exponents in many systems 0, H3, EH, 

Usui mm. 



APPENDIX B: LINEAR ENTROPY FOR 
QUADRATIC TWO-MODE BOSON SYSTEMS 

In this Appendix, we derive the expression for the lin- 
ear entropy for the ground state l^) of a general quadratic 
two- mode boson Hamiltonian of the form ^ . This linear 
entropy is defined as 



£ Vln . = 1 - Tr(p 2 A ) = 1 - Tr(pi) , 



(Bl) 



where pa,b = T^B,A (IV'XV'I)- Here, we suppose that 
bosons a and b define subsets A and B. 

As shown in Ref.pJ, the reduced density matrices pa,b 
can always be written as the exponential of a quadratic 
form in the operators related to the corresponding subset. 
Concretely, tracing out over B, one has pa = &~ K where 



K 



kq + ma) a + K2 (a) + a 2 ^j 



(B2) 



The coefficients determined from the conditions 

TrAPA = 1, ((o+o)) = (at a) and (( a t 2 )) = ( fl t 2 ), ( B 3) 

where (fl) = (ip\Q\ip) and ((0)> = Tr A (e~ K n). 
The next step is to diagonalize K which becomes 



K = No- 
where, as shown in Ref. fill ]. 

A = In- 



Ac'c, 



1 ' 



(B4) 



(B5) 



H = ^-<(at+a) 2 )((at-a) 2 ). The 



and 

condition Tt^PA = c~ E "/(l — c~ 



normalization 
1 then implies 

e --E n _ 2/(p + 1). It is then straightforward to compute 

(B6) 
(B7) 



1 _ xf e~ 2(Bo+Actc) 
1 

A 4 ' 



1 



Note that it is trivial to generalize this derivation for 
an arbitrary number of bosonic or fermionic modes using 
the expression of the reduced density matrices given in 
Ref. [l|. 



APPENDIX C: FIRST-ORDER CORRECTION OF 
THE ENTANGLEMENT ENTROPY IN A 
TWO-MODE BOSON SYSTEM 



We give the exact expression of the first-order correc- 
tion of the entanglement entropy in a two-mode boson 
system. This correction is derived in Ref. (36j for the 
general case. 

Let us consider a Hamiltonian of the following form 



H = H Q +gH 1 , 



(CI) 



where Hq is a quadratic two-mode boson Hamiltonian 
written in terms of bosons a and b, and where j « 1. 
Using the same notations as in Appendix[Bl we introduce 
the parameter 



// = y-((at + a ) ><(at-a) 
M (0) +3M (1) +CV), 



(C2) 
(C3) 



Then, as detailed in Ref. [36(, the entanglement entropy 
at order g° and g is obtained by expanding the function 



^) = ^ln^-^ln^, (C4) 



at the corresponding order. One then obtains the entan- 
glement entropy between both modes 



£ = -Tr A (pa^Pa) , 
= +g£W +0(g 2 ), 



(C5) 
(C6) 



with 

£ (0) 



^ + l, n e±l-^lln^l, ( C7) 



2 2 
— In ^ (0) ^ ] 



P 



(0) _ 1 



(C8) 



The only difficulty is thus to compute p, at order g and 
one readily has the correction of order g to the entropy. 
Note finally that the second order correction is not ob- 
tained by expanding F at order g 2 . 
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